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ABSTRACT 


A transfer function design theory for multivariable control synthesis is 
highlighted. The use of unique transfer function matrices and two simple, basic 
relationships - a synthesis equation and a design equation - are presented and 
illustrated. The basic idea of the method is straightforward, easy to understand and 
easy to apply. 

This multivariable transfer function approach provides the designer with a 
capability to specify directly desired dynamic relationships between command 
variables and controlled or response variables. At the same time, insight and 
influence over response, simplifications and internal stability is afforded by the 
method. A general, comprehensive multivariable synthesis capability is indicated 
including nonminimum phase and unstable plants. Gas turbine engine examples are 
used to illustrate the ideas and method. 


INTRODUCTION 


The concept of transfer function has been a mainstay of control engineering 
and design. A primary motivation to use transfer functions for multivariable design 
is to increase insight, choices and simplifications so that a designer may interact 
with multivariable systems in a more direct and integrated manner. The multi- 
variable transfer function idea was applied to jet engine control by Boksenbom and 
Hood (1) and Feder and Hood (2) about mid- century. Practical computation with 
transfer function matrices was a difficult issue at that time; and, the question of 
how to extend performance specifications to matrices of transfer functions posed 
another difficulty which is still not completely resolved. In this paper, the use of 
transfer functions for design of controller dynamics for linear multivariable models 
of turbine engines is reexamined. 

Control system design methods in the frequency domain traditionally have 
been of two types. In the first type, the designer works indirectly to adjust open 
loop characteristics so that, when the loop is closed, an acceptable system results. 
In the second type, the designer works directly from the closed loop specifications 
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to the specific controller dynamics required. We refer here to this second type 
method as a synthesis method. A classical discussion of both methodologies may be 
found in Truxal (3) . 

Two multivariable transfer function system relationships are derived- a 
design equation and a synthesis equation. The synthesis equation is used to display 
internally stable closed loop response possibilities; the design equation is used to 
compute and simplify explicit controller dynamics. Gas turbine engine examples 
illustrate ideas and methodology. 


MULTIVARIABLE CONTROL SYNTHESIS 


The basic notion of multivariable control synthesis with transfer functions is 
straightforward and easy to understand. Consider Figure 1, a block diagram for a 
unity negative feedback multivariable structure with no disturbances. References, 
error, plant input and plant output are designated r, e, u and y respectively. 

Assume the plant has equal numbers of inputs and outputs, thus P{s) is a square 
matrix of transfer functions. This assumption is not nearly as restrictive as one 
might suppose at the outset. More on this later. The controller G(s) is also square. 



Figure 1. Unity Feedback Structure 


The problem is, given plant P(s), to design a controller G(s) to achieve 
desired, internally stable, closed loop response T(s) as indicated in Figure 2. The 
objective is to design G(s) so that closed loop response T(s) is achieved in such a 
way that designer choices, insight and influence are made available and remain 
accessible. References (4) and (5) can provide more details for the interested 
reader. 



Figure 2. Desired Response 
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A Design Equation 


From Figure 1, the total response of the unity feedback loop is 

y = PG ( I + PG ) _1 r (1) 

The desired response is 

y = T r (2) 

Combining equations (1) and (2) and solving for G(s) gives the controller 

G = P _1 T ( I - T ) _1 (3) 

This equation may be written in a convenient, compact form 

G = P _1 Q (4) 

where Q, a performance matrix, is defined by Q = T ( I - T ) \ Equation (4) is 
named the design equation for controllers G(s) under the unity feedback structure of 
Figure 1. The design equation simply and clearly indicates that controller design 
focuses upon the properties of the plant inverse, P(s) , and how they interact with 
Q(s), the performance matrix. 

For a unity feedback structure as in Figure 1, and the case of decoupled 
response forms where the response matrix T is diagonal, a Q-T transform table 
conventently exhibits elements of the performance matrix, q^, corresponding to 
given elements of the response matrix, tjj. Table I lists some standard response 
forms and related performance element forms for unity feedback loop structures. 

TABLE 1 


A Q vs. T Transform Table - Unity Feedback, Decoupled Response 


DESIRED tji (s) 

<Ut (3) 

1 

1 

T s + 1 

T s 

1 

1/ <T X + T 2 ) 

(Tx s + 1) ( T z s + 1 ) 

T 1 T 2 

( - — ) 3 + 1 

Ti +T 2 

S: 

3 

to 

K : K = Wn /!> f T - 1 /« r«! 

s 2 + 2 f W n s + W n 2 

s (T s + 1 ) s ”n 

I 

1/(T! + T 2 + T 3 ) 

( T l S + 1 ) ( T 2 s + 1 ) (T 3 3 + 1 ) 
T 1 T 2 T 3 

3 ‘'TFt, 1 s + 1 ) ( t 3 s + 1 ) 

Ti T 2 a 3 
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More General Feedback Structure 

To effect control of a plant it is necessary to use actuators to drive inputs 
and to use sensors to measure outputs. Moreover, sensors and actuators can 
introduce significant dynamical effects into signal paths of the loop. Therefore, a 
more general feedback structure, which accommodates these effects, is shown in 
Figure 3. 



Figure 3 General Feedback Loop 

A(s) and S(s) are diagonal actuator and sensor matrices, respectively. The output 
y is sensed and becomes y s ; the input request, ur , commands the actuators to 
produce the plant input, u. 


From Figure 3, overall response of the loop is 
y = ( I + PAGHS)" 1 PAG r 
The desired response is 


(5) 


y - T r 

Combining equations (5) and (6) and solving for the controller G 

G = A -1 P -1 T ( I - HST ) _1 
A performance matrix 

Q = T ( I - HST ) _1 

is easily identified. The design equation becomes 

G = A" 1 P _1 Q 


( 6 ) 

(7) 

( 8 ) 


( 9 ) 
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This is a design equation for feedback systems depicted in Figure 3 with the 
performance matrix Q per equation (8) . Equation (9) is a generalization of equation 
(4); controller dynamics are determined by the characteristics of the actuated plant 
inverse, (PA) - *, and the performance matrix, Q. The plant inverse transfer 
function matrix is a key element in the design equation. What about the existence of 
the plant inverse? Is this a serious restriction to transfer function design? 


The Plant Inverse 

Fortunately, the need for existence of the plant inverse turns out to be not a 
significant limitation of the design equation. Rather, we can indicate to the contrary 
that the plant inverse establishes and displays vital plant characteristics needed to 
effect successful closed loop control design. Four system and plant features, 
essential for design, are established and identified by the plant inverse transfer 
function: 

1. meaningful multivariable control (6) 

2. plant trackability (8) 

3. multivariable plant zeros (9) 

4. cancellations and simplifications 

Rekasius (6) and Wonham and Morse (7) have shown that if the number of 
plant inputs equals the numbers of its outputs and if P(s) is full rank, i. e. , P(s)-1 
exists, then one has both a meaningful multivariable control problem and necessary 
and sufficient conditions for existence of a physically realizeable controller that 
decouples the system. If the number of plant outputs is greater than the number of 
its inputs, some of the outputs cannot be controlled independently. Even if the 
number of plant inputs exceeds the number of its outputs, independent control of all 
outputs is not possible if the plant transfer function matrix, P(s) is not full rank (6). 

R. J. Leake, et al (8) define a step trackable linear multivariable plant as 
one which can asymptotically achieve any constant steady-state output with a bounded 
control. It is shown that step trackability for proper rational continuous square 
plants is equivalent to the conditions that: 

1. the plant is invertible 

2. the plant has no multivariable zeros at the origin ( s = 0 ) 

Leake goes on to show the significance of step trackability by demonstrating that 
internally stable, decoupled closed loop design is possible if and only if the plant is 
step trackable. 

Importantly, the multivariable zeros of a plant, P(s), are the poles of the 
inverse, P(s)~l, Wyman and Sain (9). Therefore, multivariable plant zeros are 
readily identified from the factored form of the inverse transfer function matrix. 
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Thus, existence of the plant inverse assures conditions needed to effect 
design. Moreever, the plant inverse matrix provides essential design information 
about plant trackability, about plant multivariable zeros and about existence of 
meaningful and internally stable closed loop control realizations. 


A Synthesis Equation 

Little has been said about loop stability. From another view of the problem 
of synthesis of closed loop controllers, a synthesis equation will be derived which 
can be used to establish existence of internally stable closed loop controllers. The 
synthesis equation was first proposed by Dr. R. J. Leake of Fresno State. Connec- 
tion of the equation to current system theory (10) and internal stability was made by 
Dr. M. K. Sain of Notre Dame. The author is happy to acknowledge continuing 
collaborations and discussions with Professors Sain and Leake on multivariable 
synthesis with transfer functions. 

Consider Figure 4 where r denotes request, u denotes control action, and 
y denotes response. Under broad assumptions, there exist linear operators T: R— »Y 



Figure 4 A General Control System 

and M: R— *U, where R, U, and Y may be understood as R(s)-vector spaces of finite 
dimension such that (5) 

y=Tr, u = M r (10) 

The plant can be understood in terms of an operator P: U-+Y, such that 

y = P u (11) 

Combining equations (10) and (11) obtains the relationship 

T = P M (12) 
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Bengtsson(lO) proves that internally stable feedback realizations of systems depicted 
by Figure 4 exist if any only if M is proper and stable and T is proper and stable. 

Imposing a tracking requirement, as for example that y should asymptotically 
track step responses, then P must be epic. If in addition the number of inputs 
equals the number of outputs, then P is monic also. Therefore, we can address an 
inverse total synthesis problem ( ITSP ) (5) which is governed by the synthesis 
equation 

M = P _1 T (13) 

Note that equation (13) is similar in form to design equation (4) or (9). 


Two Basic Equations 

Two equations form a basis for multivariable synthesis with transfer functions: 

• the synthesis equation M = P”^T 

• the design equation G = P 1 Q 

The idea is, for given plant P, to select proper and stable T so that M is also 
proper and stable. This insures existence of internally stable controllers. Thus, 
the synthesis equation displays all possible responses T which have internally 
stable feedback realizations. 

Feedback realizations of M, by controller dynamics G and H, as indicated in 
Figures 1 and 3, are obtained by applying the controller dynamics design equation 
G = P -1 Q. The response matrix T maps to the performance matrix Q = PG in 
Figure 1 and Q = PAG in Figure 3. The issue of internal stability of closed loop 
realizations is still under study (4) , (5) . However, applications of the synthesis and 
design equations to numerous examples, including nonminimum phase plants, suggest 
the conjecture that if M and T are proper and stable, and if no cancellations of 
right hand plane poles and zeros occur in the open loop matrix products PAGHS 
(Figure 3 ), then internal stability of the closed loop is assured. In any event, in the 
abscence of a complete general theory and proof, doubts on internal stability can be 
resolved in practise by computer simulations of specific closed loop realizations. 

The foregoing ideas and use of the synthesis and design equations are illus- 
trated by examples. 


DESIGN EXAMPLES 


Two turbine engine examples are given to demonstrate linear multivariable 
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synthesis with transfer functions. The first design example uses a third order 
research model of General Electric's J-85 engine with two inputs and two outputs. 
The second example uses a sixth order model of Pratt & Whitney's F100 engine (11) 
with four inputs and four outputs. It is a pleasure to acknowledge the computational 
support of E. J. Olbey, S. A. Stopher and J. W. Wildrick of the Bendix Energy 
Controls Division. 


Example 1 

A transfer function matrix of the J-85 engine at sea level, 100% speed 
condition is given by 

y = P(s) u 

where y is the output vector and u the input vector. The output vector y'=(N, T) 
where N is rotor speed, RPM, and T is turbine temperature, °F. The input 
vector u' = (Wf , Aj ) where Wf is fuel flow, pounds per hour, PPH, and Aj is 
nozzle area, in^. The plant transfer function is 


5j_6 

(. 61s+l ) ( . 016s+l) 


55.7 

(. 61s+l) ( . 06s+l ) 


P(s) = 


. 17 (1.3S+1 ) 

(.61s+l)(. 23s+l) (. 016s+l ) 


-1.9 

(. 61s+l)(. 23s+l) (. 06s+l) 


Response Specification 

A closed loop controller is desired to control the engine so that system 
response to a step: 1. settles in one second, 2. has no overshoot, 3. obtains zero 
steady state error. Also, decoupling of the output is desired thus the response 
matrix T is diagonal. The control problem is pictured by the diagram in Figure 5. 
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Plant Inverse 


The plant inverse matrix is calculated first. The plant inverse exists, therefore, a 



.094 (. 016s+l ) 

.0084 (. 06s+l )(1. 3s+l ) 


2.8 (. 016s+l ) (.23 s + 1 ) 
-. 28(.06s+l )( . 23s+l ) 


meaningful problem is posed; and , it is possible to shape the response of the outputs 
independently with the available inputs. Decoupled response is possible. The plant 
has no multivariable zeroes since P - has no poles. Thus P and P - ^ indicate that 
at the given condition, the J-85 engine is a stable plant with no multivariable zeros. 
No possibility for rhp cancellations from the plant. 


Synthesis Equation 

The synthesis equation, M = P“^T, is applied to determine possible loop 
responses. Internally stable realizations exist if and only if both M and T are 
proper and stable. For diagonal T 


M = 


0.94 (. 016s+l ) 


L 11 


2. 8(. 016s+l)(. 23s+l) t 22 


M 


. 0084(. 06s+l)(l. 3s+l) t ix ~.28(.06s+l)(.23s+l) t 22 

XJL1 form. 


is proper and stable when t, ,, too = p- \ / 

F F ^ (T x s+1 ) ( T2S+1 ) (* 


Selection of Response Matrix 

To meet response specifications, the response matrix T is selected and 
structured initially as follows: 

• diagonal - decoupled response 

• predominant time constant = . 25 sec - one second settling 

• gain = 1 - zero steady state error 

• tjj = 1/ ( . 25 s ± 1 ) (*) - to satisfy synthesis equation 

The above structure of T implies that the performance matrix Q=T(I-T) 

is also diagonal and q^ = k/s (*). Of course t^ and t £2 can he chosen differently 
and independently. 
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Design Equation 


Controller dynamics are computed by the design equation G = P 
diagonal Q 



For 


.094 (.016s+l ) q u 

. 0084(. 06s+l) (1. 3s+l) q n 


2.8 (. 016s+l)(. 23s+l) q 22 

-. 28(. 06s+l)(. 23s+l) q 22 


The above matrix form clearly shows that G is simplified if q-Q = q 2 2 = K/s(.06s+l) 
(. 016s+l). Combining the q.. and t^. requirements gives 


*11 *22 


(. 25s+l )(. 08s+l)( ,016s+l ) 


and 


2. 89 


q il q 22 s ( .06s+l ) ( .0168+1 ) 
Then controller dynamics are 


G(s) = 


.27 

s (. 06 s + 1 ) 

. 02 (1.3 s + 1 ) 
s ( .016 s + 1 ) 


8. 1 (.23 s + 1 ) 
s (.06 s + 1 ) 

81 ( - 23s + 1 ) 
s ( .016 s + 1 ) 


Verification 

Computer simulation of the J-85 closed loop system (Figure 5 ) with the above 
controller dynamics verifies the desired response T = l/(.25s+l )(. 08s+l)(. 016s+l ). 
Figure 6 shows the response of the J-85 engine- control system to a 500 RPM step 
in speed request Nr only. Response specifications and decoupling are achieved. 
Figure 7 shows system response for 500 RPM step in speed request Nr and -50 
degree step in temperature request Tr. 


Example 2. F100 Turbofan Engine 

An extensive set of linear state descriptions of the F100 turbofan engine 
were given by Miller and Hackney (11) . In this example a reduced model at sea 
level, 67 degree power lever condition is controlled. This example illustrates a 
realistic design situation including engine, actuators and sensors. Use of approxi- 
mate cancellations to simplify the controller is also illustrated. 
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Engine Dynamics 


A reduced state model of the F100 engine at sea level, 67 degree power 
lever condition is 


x = A x + B u 
y = C x + D u 

where 



-4.064 

3.895 

-470.5 

7.971 

5.294 

-3.005 


.03718 

-2.958 

-59. 13 

.1727 

2.08 

12.43 

A = 

.03389 

.0067 

-4.442 

.0059 

. 1474 

.0985 

A — 

1. 164 

-2.646 

-331.6 

-50.05 

-.473 

-11.36 


.05174 

-.1176 

-14.74 

-2.001 

-2.021 

-.505 


.00184 

.0036 

-.601 

.00008 

.0009 

-.666 


.8686 

-14.51 

-96. 14 

9.246 




.9096 

-58.46 

-1.053 

-60. 15 



B a 

-.007994 

-79.66 

1.2 

.3673 



5.643 

-112.2 

-18.23 • 

41.53 




.2508 

-4.99 

-.8106 

1.846 




.01 

-.3166 

-. 02915 

.07426 





■“ 







— 

— 


l 

0 

0 

0 

0 

0 


0 ... 

0 

r* — 

0 

1 

0 

0 

0 

0 




C - 

0 

0 

1 

0 

0 

0 

D = 

• • 

• * 

• 

• 


0 

0 

0 

1 

1 

0 


• • 

• 









0 

0 


The states, x, inputs, u, and outputs, y, are 


Xj = Nj , fan speed, RPM 

x 2 = N 2 , compressor speed, RPM 

x 3 = P 7 , augmentor pressure, PSI 

* 4 = T h ^> fan turbine temperature (fast), °F 

* 5 = T i c t fan turbine temperature (slow), "F 

x 6 = T, burner temperature (slow), "F 

Uj = WF, fuel flow, PPH 


2 

u 2 = AJ, exhaust nozzle area, FT 
Ug = CIVV, inlet vane position, DEG 
u 4 = RCVV, compressor vane position, DEG 
, fan speed, nPH 
y 2 = N 2 , compressor speed, RPM 
y 3 = P 7 , augmentor pressure, PSI 
y 4 = FTIT , fan turbine inlet temperature, °F 
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The transfer function matrix of the engine (y = Pu)isP = C( sI-A ) -1 B + D and is 
shown in Figure 8. In the figure the notation (T s+1 ) = (T) and ( A s2 + b s +1 )=(a,b) 
is used to save space. A stable, sixth order plant is indicated. The reader will 
note that the poles corresponding to the time constants ( 1.43 ) and (.491 ) may be 
eliminated by approximate cancellations; thus, the essential dynamics are fourth 
order. 



. 400(1. 45)(.501)(.301)(. 132)(.009) 

. 32G(1. 24) (. 480)(. 03, . 2G6)(. 010) 

. 0025(1. 32) (. 483)(. 224)(-. 131)(021) 
. 095(1. 56) (. 470) (. 377) (. 028, . 254) 


1220(1. 43)(.404)(.302)(.0175)(-. 0004) 
190(1. 19)(,48G)(.2S3)(.U19)(-.013> 

-7. 80(1. 47)(. 501)(. 34l)(. 25G)(. 02) 

72. 6(1. 56) (. 470) (. 320)(. 196)(-. 004) 


-30. 6(1. 43)(. 493)(. 302)(, 094)(. 019) 
- 1. 90(. 027) (. 290, 1.08) (,02l)(. 012) 
. 0240(1. 95)(. 099)(. 475)(.329)(. 02) 
-1. 23(1. 55)(. 470) (. 088, . 566) (. 031) 


-10. 8(1. 41)(. 4GG)(. 300)(-. 039)(. 014) 
-19. 1(1.50)(. 433)(.027,.732)(.020) 
-.0160(1.G3){-1.4'1) ( .51) (. 122)(. 02) 
1. 83(1, 54)(.470)(. 118)(.035, . 144) 


(1.43)(.491X.302)(.02)(. 0297, .230) 


Figure 8. F100 Transfer Function Matrix 


The plant inverse matrix is shown in Figure 9. The factored form indicates that the 
inverse has two poles associated with time constants (1.55 ) and (.470 ) which are 
zeros of the plant. Again, these factors approximately cancel from the plant inverse 
matrix. 



331(1. 55)(.470)(. 109) 
00021^5(1. 54)(. 655)(. 470) 
04G(1. 55)(. 470)(. 235) 
-.0035p(l. S0)(. 470) (. 061) 


. 697(1. 5S)(. 470)(. 119) 

. 00059^(1. 54)(. 470)(-. 172) 
. 0450(1. 55)(.470)(-. 010) 

-. 0305(1. 57)(.470)(. 403) 


42. 0(1. 55) (. 470)(-. 007) 
-.129(’.55)(.470)(.097) 
-4. 50(1. 55)(. 470)(. 0003) 
-. 138(1. 50)(.470)(-. 050) 


7. 75(1. 54) (. 000)(. 019) 

. 00270(1. 41)(.4G2)(. 008) 
. 100(1. 50)(.489)(. 009) 

. 144(1. 27)(.491)(. 019) 


(1.55) ( .470) 

Figure 9. F100 Inverse Matrix 


Response Specification 

Assume the output response specifications of the F100 engine are 
1. decoupled system, 2. step response settles in 1 second, 3. no overshoot, 4. 
zero steady state error. The desired feedback structure is shown in Figure 10. 
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Figure 10. General Feedback Structure 


Actuator dynamics are given by u = A Ur and sensor dynamics are given by y s =Sy 
where 


1 

. 05s+l 


0 


0 


0 


1 

•02s+l 


0 


0 


0 


A = 


0 

0 


1 

. 2s+l 


0 


0 


1 

.ls+1 


0 

0 


S = 


0 

0 


1 

. 02s+l 


0 


0 


1 

. 02s+l 


0 

0 


0 


0 


0 


1 

. ls+1 


0 


0 


0 


1 

. 5s+l 


Synthesis Equation 

The synthesis equation M = becomes 


33(. 11) (. 05) 

hi 

. 90(. 12)(. 05) 

*22 

42(-.01)(. 05) 

*33 

7. 8(. 02) (. 05) 

*44 

-,0002(.66)(.2) 

hi 

. 0005(~. 17) (.2) 

*22 

-. 13(. 1)(. 2) 

*33 

. 003(. 01) (. 2) 

*44 

046(. 24)(. 1) 

hi 

. 045(-. 01) (. 1) 

*22 

-4. 5(. 1) 

*33 

. 16(.01)(. 1) 

*44 

-.0035(.061)(. 1) 

hi 

037(. 40)(. 1) 

*22 

-.14(-.06)(. 1) 

*33 

. 14(. 02) {. 1) 

*44 _ 


M is stable and proper when the tjj form is K/(Ts+l)(Ts+l). Based on the response 
specifications, t^ = l/(. 25s+l)(. 01s+l ) is chosen. 


Design Equation 

The design equation G = A^P -1 Q defines controller dynamics where 
Q = T ( I - HST ) _1 . Using the above T, A, S and P" 1 matrices and choosing H = I, 
the controller G(s), simplified by cancellations and approximations, turns out to be 
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G 


fK0 


-1. 2(. 11)(. 05) 
-. 001(. C6)(. 2) 
1G(.24)(. 1) 
-.013(.06)(. 1) 

u 


3.2(. 12)(. 05) 

. 002(-. 17){,2) 
. 16(-.01)(. 1) 
13(. 40) (. 1) 


152(. 05) 

4G(. 1)(. 2) 
-16 (.1) 

49(-. 06)(. 1) 


10(.02)(.05)(.5)/(. 1C) 

. 004(, 01) (. 2) (. 5)/(. 16) 
.21(. 01)(. l)(.5)/(. 1G) 

• 19(. 02) (. l)(.5)/(. 16) 


a (.01) 


Verification 

System output response of the F100 engine using the above controller was 
verified by CSMP simulations. Command responses and decoupling for a step 
request of 4 PSI Py augmentor pressure and for a step request of 50 degrees FT1T 
temperature are shown in Figures 11 and 12 respectively. 


SUMMARY REMARKS 


Linear multivariable control synthesis with transfer functions appears to be 
feasible and practical. An output response synthesis method was described using 
two basic equations both featuring the inverse of the plant transfer function matrix. 

The plant inverse matrix is key to multivariable transfer function synthesis. 
Its existence assures possibilities for plant trackability and decoupling; and, in 
factored form, it indicates plant zeros, cancellations and potential performance 
tradeoff to simplify the controller. 

Transfer function synthesis builds on classical transfer function concepts, 
is easy to understand and contacts modern theory. Features include direct design of 
output response, cancellation and approximation and insight on response adjustments 
to simplify controller dynamics. The possibility to include both sensitivity specifica- 
tions and response specifications looks promising and is under study. 

Transfer function synthesis is applicable to gas turbine propulsion system 

design. 


125 




TIME 

Figure 11. 4 PSI P7 Step 



TIME 


Figure 12. 50 Degree FTIT Step 
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